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Abstract 

We present a determination of the strange quark mass for two flavours 
(rif = 2) of light dynamical quarks using the axial Ward identity. The cal- 
culations are performed on the lattice using 0(a) improved Wilson fermions 
and include a fully non-perturbative determination of the renormalisation 
constant. In the continuum limit we find mf s (2 GeV) = 111(6)(4)(6) MeV, 
using the force scale ro = 0.467 fm, where the first error is statistical, the 
second and third are systematic due to the fit and scale uncertainties re- 
spectively. Results are also presented for the light quark mass and the 
chiral condensate. The corresponding results are also given for ro = 0.5 fm. 
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1 Introduction 



Lattice methods allow, in principle, an 'ab initio' calculation of the fundamental 
parameters of QCD, among them the quark masses. Quarks are not asymptotic 
states of QCD and so quark masses need to be defined by a renormalisation 
procedure, 

m s q {M) = Zi(M)m b -, (1) 

by giving the scheme S and scale M. 

To convert the lattice results to continuum numbers, one needs control over 
the discretisation errors and the matching relations between the lattice scheme 
and the continuum renormalisation scheme S. Discretisation errors can be kept 
small and manageable by employing an improved fermion action. But, still, the 
lattice numbers may show considerable cut-off dependence at present couplings, 
which requires that the calculations are done over a range of sufficiently small 
lattice spacings a, discretisation errors then being removed by an extrapolation 
to a = 0. The perturbative relations between renormalised quantities in the 
continuum and the bare lattice results are in almost all cases known to one-loop 
order only. Data show (for example pQ) that 0(af) corrections can be large, 
O(10 — 20%) at spacings a ph 0.1 fm, which makes a non-perturbative calculation 
of the renormalisation constants indispensible. 

Many calculations of the strange quark mass, both with rif = 2 12112111! an d 
n/ = 2+1 OiniHIE] flavours of sea quarks, employed perturbative renormalisation 
of the bare quark mass and were restricted to lattice spacings a > 0.1 fm. (Recent 
rif = 2 + 1 results [Zj, have finer lattice spacings.) These authors quote strange 
quark masses of O(80) MeV, lying 15 — 30% below the corresponding quenched 
results PHUl. 

In JTT] we have presented an entirely non-perturbative (NP) calculation of 
the light quark masses based on the vector Ward identity (VWI), using non- 
perturbatively O(a) improved Wilson fermions with rif = 2 flavours of dynamical 
quarks. The calculation was done at four different lattice spacings 0.065 < a < 
0.09 fm, which allowed us to perform a continuum extrapolation. We found a 
strange quark mass of m" s (2GeV) = 119(5)(8) MeV. A highlight and essential 
ingredient of the calculation was that we were able to compute the flavour singlet 
mass renormalisation constant, which is needed in the VWI approach. 

This result has been complemented recently by further studies in [T2*l ITS] 
which also used NP determinations of the renormalisation constant. 

In this paper we present an independent calculation of the strange quark mass 
using the axial vector Ward identity (AWI), again for rif = 2 flavours of improved 
Wilson fermions. The AWI involves only non-singlet quantities and thus provides 
an important test of our previous calculation. 

The paper is organised as follows. As we shall be considering not only the MS 
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scheme, but also the RI'-MOM scheme (which is more convenient for a lattice 
calculation) we first discuss in section |2 renormalisation group invariants (RGIs), 
taking the quark mass as an example, and how to convert to them. We also 
collect together relevant formulae for the MS scheme. Also discussed is the unit 
and scale we shall use - the ro-force scale - and thus the relevant conversion 
factor to physical units. In section El we compile some results from leading order 
(LO) and next to leading order (NLO) chiral perturbation theory (xPT) and 
re-write them in a form suitable for our calculation. Section 0] describes some 
lattice details relevant for 0(a) improved fermions. This is followed in sectional 
by the non-perturbative computation of the renormalisation constant. Relevant 
results for the RI'-MOM scheme are given, both for the the lattice computation 
of Z^~ MOM and for the conversion to the RGI form. The section is concluded 
with a comparison of this result with the results obtained by other approaches 
(principally the tadpole improved (TI) perturbation theory method). In sectional 
results are given for the strange quark mass, first at finite lattice spacing, and 
then the continuum extrapolation is performed to give our final answer. Finally 
in the last section, section [?l we compare our AWI result with the previously 
obtained VWI result and also with other recent mass determinations. In the 
appendix, tables of our raw data results for the quark mass are given. 



2 Renormalisation Group Invariants 

The 'running' of the renormalised quark mass as the scale M is changed is con- 
trolled by the (3 and 7 functions in the renormalisation group equation, defined 
by 



(3 s (g s (M) 



dg s (M) 



1 in 



(AM) 



dhgM 
dlogZg(M) 



d\ogM 



(2) 
(3) 



where the bare parameters are held constant. These functions are given pertur- 
batively as power series expansions in the coupling constant, 



-b g 3 -b ig 5 -b s 2 g 7 -bfg 9 - 
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(4) 



The first two coefficients of the /3-function and first coefficient of the 7 m function 
are scheme independent, 
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while all others depend on the scheme chosen. 

We may immediately integrate eq. ((21) to obtain 
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The renormalisation group invariant quark mass 1 is defined from the renormalised 
quark mass as 
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= AZ s m {M)m s {M) = AZ s m {M)Z s m {M)m b - 



rjBGI bare 



(8) 



where 



[AZi(M)}' 1 = \2b g s (M) 



r g s {M) 



exp 



7m (0 , ^ruO 

P S (0 bo£ 



(9) 



and so the integration constant upon integrating eq. (0) is given by A 5 , and 
similarly from eq. Q the integration constant is m^ 31 . A s and m^ 31 are thus 
independent of the scale. (Note that although the functional form of AZ^(M) 
is fixed, the absolute value is not; conventions vary for its definition.) Also for a 
scheme change S — > S' (it is now sufficient to take them at the same scale) given 
by 
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S' 
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S\2 



+ ... 



(10) 



then mf 31 remains invariant, while A changes as A 5 ' = A 5 exp(ti/(2&o))- Note 



also that analytic expressions for the integrals in eq. 
for low orders, for example to two loops we have 



or eq. (jHJ) can be found 



AZ s m {M) = 26 (/(M)) 



2b 
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bo 
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2b b 1 



(11) 



Thus we have a convenient splitting of the problem into two parts: a number, 
m^ 37 , which involves a non-perturbative computation, and is the goal of this paper 
and, if desired, an evaluation of AZ^ which allows the running quark mass to be 
given in a renormalisation scheme S. 

In the remainder of this section we discuss the evaluation of AZ^ in the 
MS-scheme, which is conventionally used, and for which four coefficients in the 



1 Analogous definitions hold for other quantities which depend on the scheme and scale 
chosen. 
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perturbative expansion are known. For the (3 function we have [T3J EH EI] , 
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and ^5 



Cs)nJ) 



1.03693. 



243 27 

C being the Riemann 



with Cs = 1.20206. . ., (4 
zeta function. 

This scheme is a manifestly perturbative scheme and so should be used at a 
high enough scale M = fi so that perturbation theory is reliable. Computing 
[A^ fS (/i)] _1 involves first solving eq. (J7J) for g MS (as a function of /x/A MS ) and 
then evaluating eq. Q. Practically we expand the j3 and 7 functions to the 
appropriate order and then numerically evaluate the integrals. The final results 
are given in Fig. ^ 

Conventionally light quark masses are defined at a scale of \i = 2 GeV, which 
means giving a value for A MS in MeV. We set the scale here by using the 'force 
scale' ?"o, which means first changing from the A MS unit to the r$ unit. From JH] 
(see also (201), we use the value 

r A^ = 0.617(40)(21). (14) 

The ro scale in MeV still needs to be set. A popular choice is 



r = 0.5 fm = 1/(394.6 MeV) 



(15) 



which is useful when making comparisons with other results for the quark mass. 
Alternatively from a fit to the dimensionless nucleon mass r m^ using results ob- 
tained by the CP-PACS, JLQCD and QCDSF-UKQCD collaborations, following 
we found a scale of 



r = 0.467(33) fm = 1/(422.5(29.9) MeV) . 



(16) 
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Figure 1: One-, two-, three- and four-loop results for [AZ^ S (fi)] in units of A 



r 


one-loop 


two- loop 


three-loop 


four- loop 


0.5 fm 
0.467 fm 


0.735(11) 
0.745(12) 


0.682(10) 
0.690(10) 


0.698(11) 
0.707(12) 


0.700(11) 
0.711(12) 



Table 1: Values of [AZ* /s (/x)]" 1 at [i 
errors in eq. (|14j) . 



2 GeV. The errors are a reflection of the 



Similar results were obtained in [221 E2]- The 0.033 fm error in the r = 0.467 fm 
estimate is roughly equal to the difference between the two r values. As these 
different tq values give an idea of the uncertainties involved in setting the scale we 
shall derive results using both values of ro and regard this as giving an estimate 
of a possible scale systematic error. The different values of r in eqs. (fl3j). (fT^j) 
give A 1 ^ = 243(16)(8) MeV and 261(17)(9) MeV respectively 

Results forJAZ^]- 1 at \i = 2 GeV are given in Table □ At /i = 2 GeV 
we have /i/A A/s ~ 8, and it seems that already at this value we have a rapidly 
converging series in loop orders. Indeed, only going from one loop to two loops 
gives a significant change in [AZ^ S ] _1 of order 8%. From two loops to three loops 
we have about 2%. The difference between the three-loop and four-loop results 
is 0{\%). So if we are given m^ 31 , and we wish to find the quark mass in the MS 
scheme at a certain scale, we shall use the four-loop result from eq. (jHJ) as shown 
in Fig. □ 
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3 Chiral Perturbation Theory 



Chiral perturbation theory (xPT) gives the LO and NLO result 2 [231201 for nf 

degenerate sea quarks of 



4tt/o, 



Xab 



1 + Xsn f (2a 6 - a 4 ) + XAs(2a 8 - a 5 ) 

1 Xa(xs ~ Xa) In Xa - Xb(xs - Xb) In Xb 



+ 



nf Xb-Xa 



, (17) 



where the first term on the RHS is the LO term, the NLO terms being the 
remaining terms and are the low energy chiral constants (LECs) evaluated at 
the scale A x = 4irf n . With our conventions for the pion decay constant in the 
chiral limit, / , we have f v = 92.4 MeV so that A x w 1160 MeV. 
Xab are related to the quark masses by [21)] . 

X* B = Bi( l A+ ™ B) \ A, Be {Vi, V 2 , S}, (18) 

where mj is the sea quark mass and m^, i — 1,2 are the (possibly) non- 
degenerate valence quark masses. In particular we have 

Xab = 1(xa + Xb), where \a = Xaa ■ (19) 
Furthermore in eq. ()18j) . Bq is related to the chiral condensate by 

Bl = -Uqq) S - (20) 
Jo 

Apart from /„., none of the other LECs (here B , or (qq), and a4,a 5 ,a 6 ,a 8 in 
eq. (|T7)l) are well determined. Typical values are (for rif = 3), a 4 = —0.76(60), 
a 5 = 0.5(6), a 6 = —0.5(4), a$ = 0.76(40) (as compiled in [21]) giving 

2a 6 -a 4 « -0.24, 2a 8 -a 6 w -1.02 , (qq) 1T ^(2 GeV) w -(267MeV) 3 , (21) 

(the (gg) result is taken from |28j). 

m^f in eq. (|T7|) is the pseudoscalar mass (with A,Be {Vi,V 2 }) implicitly 
depending on the sea quark mass mg. Again for degenerate valence quarks, we 
write i7ip S = m p/- Note that we also numerically allow for a valence quark mass 
to be equal to the sea quark mass, so for example we can write m£,. 

We shall assume eq. ()17|) as the basic functional form for the relation between 
the quark mass and the pseudoscalar mass in the following. As expected this 



2 The NNLO result has recently been constructed in |24| . 
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equation is symmetric under an interchange of the two valence quarks. However 
it does not have the most general structure allowed by this symmetry. 

Eq. (117)1 is not very convenient for comparing with numerical results for a 
variety of reasons (see section lfi~Tl and eqs. (|57jl . (|58|)). First, the quantities we 
measure are not Xab, m p^ but only proportional to them, 

AB 

Xab = c x y AB , ^ = c m Mp S B , (22) 

where dab and M p B are new variables. Substituting these into eq. (fT7j) simply 
shifts the coefficients of the various terms (including the 1 and 1/n/ terms), while 
the structural form of this equation remains the same. Note that this includes 
the cases where c x and c m are functions of the quark mass, e.g. 

c x - cf + 4% s + cf% AB + ..., 



c 



(0) 



+ c S J 1] ys + ci B %AB + .... (23) 



Secondly, we prefer to work with a function of the pseudoscalar mass rather than 
the quark mass, so we invert eq. (|17|). This gives a functional form, which we 
choose to write as 



Vab (c h - c d (l + lnc a )\ fc c + c d (l + 21nc a )\ 

Ca + h/S + VAB 



'c«A J/a(!/s - Va) lny A ~ Vb^Vs ~ Vb) lnj/ B 



c a / J/b - Va 



(24) 



Setting A = B = V this equation reduces to the degenerate valence case (and 
finally setting V = S gives the sea quark case). These sets of equations may 
be (once) iterated to produce dab I (Mp S B ) 2 as a function of (Mp S B ) 2 and (M^) 2 . 
Relevant later will be the case of degenerate valence quarks when we have 



^2 



j^y 2 = Ca + c b (M p s s ) 2 + c c {Ml) 2 + c d ((M p s s ) 2 - 2(Mj s ) 2 ) ln(M;j 2 , (25) 



which explains our original choice of the c a , Ci (i = b, c, d) coefficients in eq. 
So, as mentioned previously, we see that determining these coefficients from 
eq. ({25} which only needs degenerate valence quark masses is sufficient to find 
the results for non-degenerate quark masses, eq. (j2H - 

Numerically we shall find that higher order terms in x?T are small, i.e. c a ^> 
\ci\M 2 s (i = b, c, d) and thus all these manipulations are justified. 

The relation between c x , c m and c a is, from the LO term 

^ = c a . (26) 



8 



At NLO we have in addition 



Cd 

C n 



(27) 



and for the remaining LECs the relationships 



2oiq — 0:4 
2a 8 - a 5 



1 



n 



f 
1 



111 c b c d 
l + ln?if h in — 

Cd c a 



1 + 2 Inn* + — + 2 In — 



c d 



(21 



We write the results only for the case when c$ = c^> = (see eq. (|23p). This is 
sufficient for the later estimation of the LECs in section El as we shall be using 
continuum results 3 . 

Finally we have to find a formula for the strange quark mass. We first note 
that we have 

• Two degenerate sea quarks m s = m ud = \(m u + 

• Three possible valence quarks, m u , rrid and m s (where V denotes the 
strange quark). We write m u = m u d — Amy and vn d = m ud + Am u d 
where Amy = (ma — m u )/2 is proportional to the difference between the 
down and up quark masses. 



Inputting meson data we use the K + (us), where we set A = u, B = s in eq. 
K° (ds), where we have A — d, B — s, and together with the 7r + (ud), with A = u, 
B = d, gives after some algebra the result 



Ml + + Ml 



M 2 K+ 



LK+ "I" 1V1 K° 

+ (c b - c d ) 

+ \(Cc + C d 

-(c b + c c )M* 



+ M 2 



M 2 K+ + Mf 
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-c d 



M 2 K+ + M 2 



K° 



M 2 K+ 



ML 

2 



+ M 2 



K" 



Ml 



In (M 2 



K-* 



+ MU - M 2 



+c d M* + lnM 2 + + ---, 
for the strange quark and 

Vud = c a Ml + + (c b + c c )M* + - c d M*+ lnM 2 + + 



(29) 



(30) 



J For the more general case in eqs. (|26[1 . (I27f) and 1281) 



are replaced by c™ and c 

(0K 2 w S(l) 



JO) 



respectively, together with additional terms on the RHS of eq. Ij28(l of — l/(n/(c^- ) 2 )(c 

2cm^ /cm) and — l/(c^) 2 (c^ B ^ 1 ^ — 2cm B ^ /cm) for the second and third equations, respec- 
tively. 
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for the light quark, where the • • • include higher order terms in (i.e. NNLO) 
and terms of 0((Am u d) 2 ). 

The results of eqs. (J2*9*j) and (pffij) are valid for 'pure' QCD. To include elec- 
tromagnetic effects, we use Dashen's theorem which says that if electromagetic 
effects are the only source of breaking of isospin symmetry (i.e. m u = m^), the 
leading electromagnetic contribution to m? K+ and m 2 + are equal while m^ and 
m 2 K o are unaffected (see e.g. [20! )• Thus the masses in eqs. (j2HJ) and (IHUjl may be 
written as [30j 

m 2 K+ = Krr-Krr+m 2 , 

m 2 K o = (m£T) 2 , 
™l + =ml = (m^f- (31) 

where m^+, m K o, m n + and m n o are the 'pure' QCD numbers, while m 1 ^ 1 ', m^, 
m 1 ^ T and m^ PT are the experimentally observed numbers. Dashen's theorem 
has corrections 0(a QED m q ) from higher order terms in xPT, estimates vary as 
to the magnitude of this correction [3HE2], see [2011331 for a discussion. Recent 
results from the lattice approach [3U [3H] seem to indicate only a mild breaking 
of Dashen's theorem. 

Note that to LO in xPT using experimental values of the 7r and K masses, 
namely 

m EXFT = Q m EXFT = Q ? 

mfT = 493.7 MeV mfT = 497.7 MeV , (32) 

we have the result 

mf (M) m^-+ + m^o - m 2 + 

= {mK+ } + ( ^ 2 ( - + } « 25.9 , (33) 



independent of the value of c a . So if we are in or close to this regime, once we 
have determined the strange quark mass, this immediately gives an estimate of 
the light quark mass. Incorporating the NLO terms needs a determination of all 
the c a and q {i = b,c,d) coefficients in eqs. (J201), (1301) an d gives the results in 
section I 



4 The Lattice Approach 

Here we shall derive results for the unquenched (rif = 2) strange quark mass 
using the axial Ward identity. All our numerical computations are done with 
degenerate valence quark masses. 
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The starting point is the AWI; in the continuum we have the renormalised 
relation 

d^ = 2m s q (M)T s (M) J (34) 

where A n and V are the renormalised (in scheme S) axial current and pseu- 
doscalar density respectively. In the (bare) lattice theory the current quark 
masses are also defined via the equivalent AWI 4 

d^A, = 2m q V + 0(a 2 ), (35) 

where A and V are the 0(a) improved unrenormalised axial current and pseu- 
doscalar density 

A„ = (l + b A am q )(A ll + c A adlf r P), 
V = (l + b P am q )P, (36) 

with 

K = ll^q , P = qisq • (37) 

Wilson-type fermions allow several different definitions of the quark mass. We 
denote the (bare) quark mass defined from the AWI with a tilde, am q , while that 
from the VWI is given by 




where n q is the Wilson hopping parameter, defining the quark mass (both sea 
and valence). The critical sea quark hopping parameter, K qc , is defined for fixed 
(3 = 6/(7q (where g is the lattice coupling) by the vanishing of the pseudoscalar 
mass 5 , i.e. m„ s L - K s = 0. 

Returning to eq. (}3l)j) . ca is known non-perturbatively [55] . but not 6a and bp. 
However results using one-loop perturbation theory [57], or for non-perturbative 
quenched QCD [38J, show that the difference b A — bp is small (there is however 
an increase between the perturbative and quenched non-perturbative results). 
Multiplying by am q thus gives a correction of perhaps half a percent, which with 
our present level of accuracy we can ignore. 

Forming lattice correlation functions means that the quark mass can be de- 
fined and determined from the ratio 6 



2(P(t)O(0)) 



4 9^ 4T is the symmetric lattice derivative, conventionally chosen to be (dff 1 ^ f)(x) — [f(x + 
afi) — f(x — aft)]/ (2a), where fi is a unit vector in the \jl direction. 

5 We shall suppress the 'V superscript on the pseudoscalar mass and only include an 'S" 
superscript where necessary. 

6 Note that to reduce noise, derivatives of operators on the lattice are taken as compact as 
possible, consistent with the given symmetry. Thus we use, no fi summation, (d^ 1 ^ f)(x) = 
[f{x + afi) - 2f(x) + f(x - af L )]/(2af = (d^d^f)(x) + 0(a). 
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Table 2: The /?, and c sw values 8 and the lattice volume V = N s x iVy. The 
collaboration that generated the configurations is given in the last column. 

„ (grA 4 (t)0(o)) ^ j^p(t)o(o)) , 

2(V(t)O(0)) A 2(V(t)O(0)) { ' 

= amf + c A am[ l) + 0(a 2 ) , (39) 

where ~ in the second equation signifies that we have dropped the correction 
factor 1 + (6^ — bp)am q . O is an operator with a non-zero overlap with the 
pseudoscalar particle. We choose it here to be p smeared ; where we have used 
Jacobi smearing (see the appendix) on the operator. 

The parameter space spanned in our numerical simulations is given in Table El 
The notation is standard for the parameters of the action, see for example [40J. 
(The critical Wilson hopping parameters, K qc , have been determined in JT] for 
each (3.) 

In the appendix we list our K q for each n q together with the corresponding 
partially quenched am ps and bare AWI quark masses am q °\ am^ 1 ' and am q . 



5 Renormalisation 

5.1 Generalities 

Imposing the AWI on the lattice, eq. (J3SJ), up to cut-off effects means that the 
axial current as well as the pseudoscalar density and the quark mass must be 

For the number of trajectories generated for each for example |39| . 
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renormalised. Thus we have 



A* = Z A A^ , V s \M) = Z S P {M)V , (40) 



giving 



or in RGI form 



m s q {M) = Zi{M)m q , Zg(M) = , (41) 



<* = Z *™« > = ^(M)Zg(M) . (42) 



5.2 Non Perturbative renormalisation 

We shall employ here the RI'-MOM scheme jlj, which is easily transcribed to 
the lattice. Our implementation of this method is described in j32] ■ As discussed 
in section |2] to obtain the RGI quark mass, we must determine both AZ™'~ MOM 
and Z™ ~ MOM . The RGI quark mass can then be easily converted back to the MS 
scheme. 



5.2.1 



/\ ym'-MOM 



We start with AZ™'~ MOM . To write down the perturbative expansion, a definition 
of the coupling constant is required. The anomalous dimension coefficients have 
been determined to fourth order in jU UHl S3! by taking the coupling constant 
to be g MS . Thus the anomalous dimension function is considered as a function of 
g MS . (Other definitions of the coupling constant are possible, more closely related 
to MOM schemes ESI.) So we write 



Rl'-MOM f „MS\ 
Jm 



d m0 (g MS ) 2 + dZ{ MOM {9 



RI'-MOM I „MS\4 



(43) 



with coefficients given by 

2 



Rl'-MOM 



JRl'-MOM 
a m2 



JRl'-MOM 
a m3 



(4vr) 4 
2 

(4tt) 6 
2 

(4tt) 8 



126- —n f 
9 

20174 3344 A 47588 128 >N 856 2 ' 

G — C^iTif H Tit 

3 3 S3 V 27 9 W ' 81 f . 

141825253 7230017 , 6160 , 



324 54 

,3519059 298241 



54 



27 



<3 



,611152 5984 > , . 
+ (-rr^ — Ca)n 



243 



27 



4160 

16024 , 
729 f 



(44) 
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Figure 2: [AZ r f- MOM i^)}- 1 versus /i p /A MS . 



which allows [AZ^f MOM {Hp)\ 1 , where \i v is taken to be the momentum scale in 
this scheme, to be computed in the usual way, 



[A^'Wr 1 = 2b Q g MS ^ p f ^ exp / 



7^1 to 

/3 MS (0 



(45) 

(This result may also be shown by changing the integration variable in eq. (J0J) from 
some defined g™'- MOM to g^, by using eq. (HUJ) , i.e. g m '- MOU {n P ) = G(^(/i p )).) 

In Fig. El we show [AZ^'^™]" 1 as a function of fx p /A MS . The convergence 
for AZ"' _MDM seems slightly worse in the region of interest than that for the 
corresponding AZ^ S as there is more of a change from the 2-loop to 3-loop result 
(see Fig. ^ for a comparison with AZ^ S ). However this is mitigated by the 
extraction of RI'-MOM being performed at a range of scales including higher 
scales than 2 GeV. Thus, for example, a typical value of (a/x p ) 2 ~ 4 (see Fig. HJ) 
corresponds to fi p /A MS > 20 where the convergence (between the 3-loop and 
4-loop results) appears to be better. 



5.2.2 Z m 



RI'-MOM 



The RI'-MOM scheme considers amputated Green's functions (practically in the 
Landau gauge) with an appropriate operator insertion, here either A or P. The 
renormalisation point is fixed at some momentum scale p 2 = fi 2 , and thus we 
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have, e.g. [HJII2], 



zg'- MOU (n P ) q - ¥L 



(46) 



where Tq are one-particle irreducible (1PI) vertex functions, and Z q is the wave- 
function renormalisation. We have generated To only at values of the sea quark 
mass, i.e. the values in Table El 

To obtain Z^- MOM we need both Z A and Zf- MOM in the chiral limit. Z A is 
unproblematical, we make a linear fit of the form 

Z A = A A + B A am q , (47) 

where am q is defined in eq. ()38|) and has been determined in [TTj . 

From the Ward identity obeyed by Tp, we expect, due to chiral symmetry 
breaking, that Zp ~ MOM develops a pole in the quark mass. Hence we have Z P — > 
with the quark mass. Thus following jlH] we try a fit of the form 9 

[Zf- MOM Y x = A P + — . (48) 

am g 

From the operator product expansion [13 EHj we expect 

Bp( "> )K (^r (49) 

where the constant of proportionality is proportional to the chiral condensate. 
Thus we see that as the scale increases, the Bp coefficient decreases. 

A chiral extrapolation using eqs. (|4T|) and (}4*%|) has been made using only the 
sea data sets to determine the functions A A (/i p ) and Ap(/i p ) respectively. The 
lattice momenta originally chosen for the heaviest quark masses were kept fixed 
for the extrapolation over the different masses. If another data set did not have a 
particular momentum a linear interpolation was performed between the adjacent 
momenta straddling the given momentum. 

Results for zf' MOM are shown for (3 = 5.20, 5.40 in Fig. B P is small 
and decreasing with increasing (a/i p ) 2 as required from eq. ()49j). (Note that 
numerically, there is evidence for a term of the form Bp/am q as attempting a 
linear extrapolation, as for Z A in eq. ([47)1. gave a substantial increase in the fit 
X 2 .) Although we do not determine the chiral condensate in this section (see 
section lfi~Tj) . we note that numerically the Goldstone pion contamination to Tp 
appears to be small [5T)] . 

9 It made little difference to the Ap (or A A ) coefficients whether afh q or (am ps ) 2 is used. 
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Figure 3: A P and B P for Z£ 



Rl'-MOM 



for = 5.20, and /3 = 5.40. 



/3 


7 RGI 


5.20 


2.270(12) 


5.25 


2.191(24) 


5.29 


2.177(14) 


5.40 


2.124(06) 



Table 3: Values of ZJP from the NP method of section l5~2~T 



5.2.3 Zr 



Zf~ MOM , forming the ratio Z A jZf~ 



Taking these values for A a = Z& and A p 
and multiplying by AZ™'~ MOM will then give Z^ 1 . This should be independent of 
(afip) 2 . Some results are shown in Fig.HJ Due to cut-off effects, non-perturbative 
contributions etc., they are not quite constant although the curves become flatter 
for increasing (3. To allow for the non-constant remnants we make a phenomeno- 
logical fit of the form 



F(an P ) = r 1 + r 2 (an P ) 2 + 



r 3 



(50) 



where we associate Z^ 1 with n. This gives the results in Table El We start the 



fit range at [a\i v 



1.5, the fit results for t\ were found to be insensitive to 



decreasing the fit range. 
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Figure 4: Z^ 1 for = 5.20 (filled circles), = 5.25 (filled squares), = 5.29 (filled 
upper triangles), = 5.40 (filled lower triangles) together with fits as in eq. (|50|) . 

5.3 Comparison of Z^ 1 with other results 

As many computations of the strange quark mass have used tadpole improved 
perturbation theory together with a boosted coupling constant for the determi- 
nation of the renormalisation constant, it is of interest to compare our results ob- 
tained in the previous section with this approach. Our variation of this method, 
tadpole- improved renormalisation-group-improved boosted perturbation theory 
or TRB-PT, is described in |51| 140 j. Here we recapitulate the method. Regard- 
ing the lattice as a 'scheme', then from eq. (jHJ) we can write 

mf = AZ*f(a)m g (a), (51) 

where the renormalisation-group-improved AZ 1 ^ tr (a) is given by eq. Further- 
more in this 'lattice' scheme, we choose to use = ^q/^oc where Uq = (^TrU n ) 
(U D being the product of links around an elementary plaquette) rather than 
go, as series expansions in g D are believed to have better convergence. This is 
boosted perturbation theory. (We shall use chirally extrapolated plaquette val- 
ues as determined in [TH] at our values and so we add a subscript 'c' to u .) 
In the tadpole-improved, or mean field approximation, renormalisation constants 
for operators with no derivatives are oc uq c , which indicates that Z^'u^ will 
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5.20 


1.837 


5.25 


1.851 


5.29 


1.862 


5.40 


1.891 



Table 4: Values of I( - TRB ^ from section 15.31 



converge faster than Z^ 1 alone so we re-write the two-loop equation eq. (fTTj) as 



10 



= A ^m (a) = «0c 26 # 



2b 



6l r, 

00 



26 i>l ' 4 iij 



(52) 



where p\ is the first coefficient in the expansion of uo c , i-e. uq c = 1 — \g^pi + ■ ■ ■ 
with pi = g. 

It remains to determine d^f. This may be found by relating the (known) 
perturbative result for Z^ s to AZ^f (for simplicity at the scale \i = 1/a) by 

Z% s (\la) = - Jl} \ = 1 - g 2 nB™ s (l) + ... , (53) 



B™ s (c sw ) = (-6.79916 + 2.4967c sw - 4.28739cM • (54) 



where 

(4tt) 5 

This result for B^ s = Ba — B f is is taken from 52 . (Indeed we could alternatively 
consider TRB-PT for Z& and Zp s separately and then form the ratio. The results 
turned out to be about 1% lower than those presented here.) 

Expanding the ratio in eq. (f53*j) . by using the results in eqs. (J52J) and (fTTj) 
for the 'lattice' and MS schemes respectively and t\, defined in eq. (fTOj) (which 
using the notation of [JJJ] is numerically given by t\ = tf 1T (l), with tf 4T (c sw ) = 
0.4682013 - 71/(0.0066960 - 0.0050467c sw + 0.0298435O) S ives finall y 

C" = <C + - Pi) - 2b B^(l) = • (55) 

Thus from eq. (fH2*jl various values of AZ^ , or equivalently Z^ ;I ^ TRB ~ FT \ can be 
found. Results are given in Table |H 

We now turn to a comparison of the results. In Fig. we plot Z^ 31 versus 
(3. Our NP results from section 15.21 are shown as filled circles. They are to 

10 The TRB-PT subscript in brackets is there only to distinguish the results from those 
obtained in section 15^1 
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Figure 5: Z^P 1 versus (3. The black circles are the results from Table |3 while the 
open squares are the TRG-PT results from Table 0J Furthermore the open diamonds 
and triangles are the NP results from ^2]> using the two different results for the axial 
renormalisation constant, [53]. (The empty triangle results has been slightly displaced 
for clarity.) 



be compared with the TRB-PT results denoted by empty squares. While there 
is a difference between the results, it is decreasing for (3 — > oo and thus may 
be primarily due to remnant 0(a 2 ) effects, which disappear in the continuum 
limit. That various determinations of Z^ 31 have different numerical values can 
be seen from the results of ^2] (open diamonds and triangles). In these results 
two different definitions of the axial renormalisation constant have been used, 
|q"3] . Za is computed when dropping certain disconnected diagrams, while Z^ 
includes them. The difference between the two definitions is an 0(a 2 ) effect. 
Using Z^ OJV in Z^ 31 leads, perhaps coincidently, to very similar results to our NP 
results. 

Investigating the possibility of 0(a 2 ) differences a little further, we note that 
if we have two definitions of Z^ 1 then if both are equally valid, forming the ratio 
should yield 

yBGl(X) 

Ri = ^ r = l + 0(a 2 ), (56) 

where Z^ 1 is the result of section 15.21 and X is some alternative definition (i.e. 
TRB-PT, ALPHA-Z4, ALPHA- Zf N ). In Fig. H we plot this ratio for these 
alternative definitions. The r /a values used for the x-axis are found by extrap- 
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Figure 6: versus {a/r ) 2 c for X = TRB-PT (open squares), X = ALPHA-Z^ 
(open diamonds) and X = ALPHA- Z^ N (open triangles). 



olating the r /a results to the chiral limit. This extrapolation and results for 
( r o/ a )c ar e given in [TUj . 

We see that (roughly) all three ratios extrapolate to 1 which implies that 
any of the four determinations of Z^ 1 may be used. This includes the TRB- 
PT result. Of course other TI determinations might not have this property, and 
also their validity always has to be checked against a NP determination, so this 
result here is of limited use; it is always essential to make a NP determination 
of the renormalisation constant. Furthermore it is also to be noted that differ- 
ent determinations can have rather different 0(a 2 ) corrections, so a continuum 
extrapolation is always necessary. 



6 Results 
6.1 Generalities 

In Figs. |Hl Eland HH we plot the ratio rQm 1 ^ 31 / (r m ps ) 2 against (r m ps ) 2 , where 
r = Tq . (i.e. r depends only on the sea quark mass, S). The numerical values 
of r /a are given in [TUj . Table 2. r seems to be a good scale to use because this 
ratio numerically does not vary much, as can be seen from the figures. 
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scribed in the text. Results for equal sea and valence quark masses are denoted by an 
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p 


r RGI 


n BGI 
C b 


r BGI 


r .BGI 
C d 


5.20 
5.25 
5.29 
5.40 


0.1115(53) 
0.1169(42) 
0.1166(29) 
0.1218(24) 


-0.00227(187) 
-0.00428(160) 
-0.00199(117) 
-0.00324(079) 


0.00557(491) 
0.00612(399) 
0.00470(287) 
0.00646(179) 


0.00121(133) 
0.00153(111) 
0.00121(090) 
0.00189(056) 


oo 


0.1330(74) 


-0.00378(254) 


0.00789(607) 


0.00275(179) 



Table 5: Values of c^ 31 and cf 31 (i = 6, c, d) together with their extrapolated (contin- 
uum) values (P = oo). 

Using eq. (j2*3|) (the case of degenerate valence quarks) we set 

y v = r m^' , = r m s ps , M^ s = r m ps , (57) 

to give a fit equation of the form 

^^- 2 =cf + c™(r m s ps ) 2 + c™(r m ps ) 2 
{r m ps y 

+ cf 1 ((r m s ps ) 2 - 2{r m ps ) 2 ) \n{r m ps ) 2 . (58) 

We use this equation to determine the coefficients c^ 31 and cf 31 , i = b,c,d. The 
coefficients depend only on lattice simulation quantities and the unit chosen, and 
not on the scale as given for example in eq. (JT5j) or eq. (|16p. This can be useful, 
as an aid, when making comparisons with other results. 

As already mentioned, up to this order in chiral perturbation theory no input 
from lattice simulations with non-degenerate valence quark masses is needed. 
To test also numerically that the effects from non-degenerate quark masses are 
indeed small, we calculated correlation functions using non-degenerate valence 
quark masses ^ as well as degenerate valence quark masses with = 
{rriq +m,q)/2. We found the relevant quantities am ps and am q to differ by < 1%. 

In Figs. 171 IHll9landfTUlthe dashed lines (labelled l ^/2m K r ) represent a fictitious 
particle composed of two strange quarks, which at LO gives from eq. (}2*9"|) 
(or equivalently eq. (jSHJ)) the line (r m ps ) 2 = (r m K +) 2 + (r m K o) 2 — (r m 7r +) 2 , 
while the dashed-dotted lines (labelled 'm^') represent a fictitious pion with mass 
degenerate u and d quarks given from eq. (jHOjl (or equivalently eq. ()64))) by 
(r m n+ ) 2 . 

The presence of a chiral logarithm in the data manifests itself in the bending 
of the results for smaller quark mass, which can be seen in Figs. El EH an d El 
Results for the fit parameters are given in Table 03 

Numerically we expect the leading order in x?T to be dominant with the NLO 
giving only minor corrections 11 , i.e. we expect c^ 31 3> {rQm ps ) 2 cf 31 (i = b, c, d) and 
this is indeed found. Using these results for c^ 31 , cf 31 (i = b,c,d) in Fig. ITT1 we 

11 Alternative plots, using constant a so arh q / \am ps ) 2 against (am ps ) 2 or equivalently 
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and Cj 31 (filled diamonds) versus (a/ro) c 2 . The (ro/a) c i.e. the chirally extrapolated 
values for r^ja are used. The open symbols represent the values of c^ 31 and cf 31 
{i = b, c, d) in the continuum limit. 



present their continuum extrapolations. The r^/a values used for the x-axis are 
found by extrapolating the previously used ro/a results to the chiral limit. This 
extrapolation and results (for (r /a) c ) are also given in JH]- This gives values in 
the last line of Table El (/? = oo). 

As a first check on these results and to find some idea of possible systematic 
effects, we have varied the fit interval from the chosen (rom ps ) 2 < 5 to (rom ps ) 2 < 
4 or 6 or oo, i.e. include all the data. (We can go no lower than (rom ps ) 2 < 4 if we 
wish to have at least two sea quark masses in the fit for each (3.) There was little 
change in the fit result and we shall take these chang systematic error, see 

below and section E21 Using (rom ps ) 2 < 5 means that our pseudoscalar masses 
range from about 440 MeV to about double that value. 

As a second check on the validity of these results, let us relate them to the 
LECs in eqs. (|26jl. (|2*7j) and (|28|) (in the continuum limit). From Table El and 
eq. (J2B1) we find that 

2a 6 - a A « -0.21 , 2a 8 - a 5 w -1.19 . (59) 

These numbers are to be compared with the results of eq. (}2*T|) namely ~ —0.24 

row,/ {romps) 2 against {rom ps ) 2 using (ro/a) c which is the chirally extrapolated ro/a would 
give larger NLO corrections. 



24 



and ~ —1.02 respectively. Furthermore we have from eqs. (fTHjl. (JTHJ), (J2~2~j) and 
(loT|) the relations 

2B™ 31 1 
Cx " (4vr/ ) 2 r ' ° m ~ 4^ ' (6 ° } 
where 5^ = [AZ^(M)]- 1 ^ , which give 

£<T = — *— • (61) 

Together with eq. (J2D1) (using / « f v ) and [AZ^f (2 GeV)]" 1 from Tabled this 
leads to (cf eq. (f2T])). 

(aa\^(?r,V\- / (263(5)(5)(5)MeV) 3 for r = 0.5fm 

Wg) (2GeVj ""\ (267(5)(5)(5)MeV) 3 for r = 0.467 fm ' {b2) 

where as discussed in section we consider two scales r$ = 0.5 fm and r$ = 
0.467 fm. The first error is statistical and the second is systematic ~ 5MeV 
determined by the change in cf 31 when changing the fit interval, as discussed 
above. The third (systematic) error is due to the choice of r$ scale, taken here 
as the difference between the results and also ~ 5 MeV. 

We find it encouraging that the cf and cf 31 (i = b, c, d) results from the fit 
give numbers in rough agreement with phenomenological expectations for 2a$ — 
«4, 2«8 — «5 and the chiral condensate (and also with other lattice determinations 
of the chiral condensate, e.g. jlUlEII). However to obtain more accurate results 
will require much more precise numerical data 12 . 



6.2 The quark masses 



We now turn to the evaluation of the strange (and light) quark masses. After 
finding the coefficients cf, cf 31 {i = b, c, d), these can be substituted into eq. (I21JI) 
to give for the strange quark mass 



r m 



+ c 



(r m K +) 2 + (r m K o) 2 - (r m n + 
(r m K +) 2 + (r m K o) 2 



+ \{cf + cf 1 ) (r m K+ )' 2 + (r m K of 



2 



-(cf + cf 



RGI 



-ru c )(r m n +) 

(r m K +) 2 + (r m K o) 2 (r m K +) 2 + (r m K o) 2 - (r m lv +) 2 
x In ((r m x +) 2 + (r m K o) 2 - (r m n +) 2 ^ 



+cf I (r m 1T +f ln(r m w +) 



(63) 



12 It is possible to obtain expressions and hence in principle results for fo, (qq), 2ag — ct4 and 
2a s — a5 in terms of cf and cf 31 {i — b, c, d) . However these then all depend on the less well 

RGI (A _ 



determined NLO cf 31 (i = b, c, d) 
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Similarly for the light quark mass we have 
rom^ = cf'irom^ f + (c£* + cf* )(r m^) 4 - cf 1 (r m„ + ) 4 Hr m n+ ) 2 . (64) 

We first consider the strange quark mass. As can be seen from Table |S] or 
Fig. ^2 the errors of the NLO parameters, i.e. cf 31 (i = b, c, d) are the same size 
as the signal itself and thus using them directly in eq. (|6l3j) simply gives a change 
in the LO result (i.e. using only the c^°' term) of a few percent, together with 
a similar increase in the error (especially using error propagation; note that the 
third and fifth terms in eq. (jOBI give the main contribution to the NLO term). 

To reduce the total error on the result, it proved advantageous to use eq. (JHUI) 
to eliminate c^ 31 from eq. (|58|) in terms of 

a ' (r m K +) 2 + (r m K o) 2 - (r m 7T +) 2 ' 
This results in a modified fit function of the form 



RGI 



cf 1 + cn(r m* ) 2 - d b ] + cf'^m^f - d c 



{r m ps 



((r mpj 2 - 2(r m ps ) 2 ) ln(r m ps ) 2 - d d , (66) 



where di {i = b, c, d) can be read-off from eq. (JfiHj) and have the effect of shifting 
the various terms in the fit function by a constant. For example, the simplest to 
evaluate, c4, is given by (r m 7r +) 2 . Note that the fit coefficients cf 31 {i = b,c,d) 
are unchanged from those given in Table El Also the numerical values of the fit 
function are unchanged and are given by the curves in Figs. EJ El El and EH 

Note that, although we now have to choose the scale before we make the fit, 
the advantage is that the error on cfP 7 gives directly the error on the strange 
quark mass up to NLO. Given this, it is no longer a disadvantage to consider the 
continuum extrapolation for m" s (2GeV) directly. 

As discussed earlier, we shall consider two r scales. Given [AZ^ S (2 GeV)] -1 
from Table ^ and using the experimental values of the ir and K masses, eq. (p?2*|) 
to determine the 'pure' QCD pseudoscalar masses in eq. (JIITj) . we find the results 
in Table El 

A continuum extrapolation (using r = 0.5 fm) is shown in Fig. fT2l together 
with a comparison with our previous VWI results [TTj . The two methods have 
different 0(a 2 ) discretisation errors, but should agree in the continuum. For the 
VWI method an extrapolation gave results of 126(5) (8) MeV, 119(5) (8) MeV for 
ro = 0.5 fm and 0.467 fm respectively. Consistent agreement between the AWI 
and VWI methods within error bars is found. Finally we quote our result for the 
strange quark mass 

m Ms (9Cp v\-l H7(6)(4)(6)MeV for r = 0.5fm 

m s (z^evj-j 1U ( 6 )(4)( 6 ) MeV for rQ = .467fm ' {0i> 
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p 




rof s (2GeV) 




ro = 0.5 fm 


r = 0.467 fm 


tq = 0.5 fm 


r = 0.467 fm 


5.20 
5.25 
5.29 
5.40 


0.1179(21) 
0.1233(20) 
0.1216(20) 
0.1282(18) 


0.1175(21) 
0.1231(20) 
0.1214(21) 
0.1281(18) 


98.71(2.33) 
103.27(2.33) 
101.79(2.33) 
107.32(2.26) 


93.35(2.29) 
97.77(2.30) 
96.41(2.31) 
101.73(2.25) 


oo 


0.1393(46) 


0.1395(46) 


116.5(5.6) 


110.7(5.5) 



Table 6: Values of c^ 37 and m^ ts (2GeV) together with their extrapolated (continuum) 
values (P = oo). cf 31 (i = b, c, d) are given in Table |SJ 



where the first error is statistical and the second is systematic pa 4MeV. As 
discussed in section I6~T| we have determined it from the effect on c^ 37 of changing 
the fit interval. Furthermore the additional third (systematic) error due to the 
r scale uncertainty is pa 6 MeV, see section |21 

As discussed earlier, from the values of the c^ 31 , cf 31 (i = b, c, d) coeffi- 
cients we know that the corrections due to the NLO terms are small. Using 
the continuum value of c^ 31 from Table 03 gives for the LO strange quark mass 
mp(2GeV) = 111(6) MeV and 106(6) MeV for r = 0.5 fm, r = 0.467fm re- 
spectively. The difference is about 6 MeV, which means that NLO terms give 
about a 5% correction. We have not tried here to estimate the effects of higher 
order terms in xPT, [24J. 

For the light quark mass, the numerical situation is more fortunate. From 
Table EJ we see that |((r m 7r +) 2 (cf ;;/ + cf 37 ))/^ 1 ] ps 0.004 and similarly for 
|((r m 7r +) 2 ln(r m 7r +) 2 c I f ;/ /cf G ' I )| ps 0.005, so corrections from LO to NLO %PT 
are at the 1% level and are negligible here. We shall just quote the LO result of 

m —s (9CpV ,_\ 4.30(25)(19)(23)MeV for r = 0.5fm 
ud{ > \ 4.08(23)(19)(23)MeV for r = 0.467 fm ' K > 

where again the second error is systematic. The third (systematic) error is due 
to the scale pa 0.23 MeV. 

Finally, because the NLO corrections of x?T are small, we see that the ratio 

mf(2 GeV) 



<f(2GeV) 
is close to the LO result, eq. (J3HJ) . 



27.2(3.2) , (69) 



7 Comparisons and Conclusions 

In this article we have estimated the strange quark mass for two flavour QCD and 
found the result in eq. ([67)1 . using 0(a) improved clover fermions and taking into 
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Figure 12: m" s (2GeV) versus (a/ro) c 2 (filled circles), together with a linear contin- 
uum extrapolation (empty circle). The scale used is ro = 0.5 fm. For comparison we 
also show the results from using the VWI, open squares, [TT] . 



consideration non-perturbative (NP) renormalisation, the continuum extrapola- 
tion of the lattice results and the use of chiral perturbation theory. The NLO 
chiral perturbation theory yields a correction of about 5% to the LO result, and 
the relevant low energy constants are in rough agreement with the phenomeno- 
logical values. 

It is also useful to compare our results with the results from other groups. 
In Fig. we show some results for rif = 2 and n/ = 2 + 1 flavours (keep- 
ing the aspect ratio approximately the same as in Fig. H2*|) . A variety of ac- 
tions, renormalisations, units and scales have been used (so the results have been 
plotted in physical units using the authors' preferred values). In particular the 
HPQCD-MILC-UKQCD and HPQCD [U] collaborations use improved stag- 
gered fermions. These fermions having a (remnant) chiral symmetry and are in 
the same situation as overlap/domain wall fermions where there is no distinction 
between VWI and AWI quark masses; the bare quark mass in the Lagrangian 
simply needs to be renormalised. 

As seen earlier in section 1531 it is noticeable that the (tadpole improved) per- 
turbative results lie lower than the non-perturbatively renormalised results. Also 
results with a < 0.09 fm (i.e. a 2 < 0.008 fm 2 ) appear to be reasonably consistent 
with each other (this is more pronounced for the AWI results than for the VWI 
results). While results for a < 0.09 show some lattice discretisation effects, using 
results at larger lattice spacings seems to give a fairly constant extrapolation to 
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Figure 13: Results for m" s (2GeV) versus a 2 using the AWI (upper plot) and VWI 
(lower plot) methods. The results are presented with the collaborations preferred units 
and scales. Circles (together with a linear continuum extrapolation) are from this work 
and diamonds from squares from ^3]; up triangles from [3]; down triangles 
from [1]; left triangles from jSJ; right triangles from [S]. NPR denotes non-perturbative 
renormalisation, while TB-PT denotes TI boosted perturbation theory. [HI IB] are for 
rij = 2 + 1 flavours; the other results are all for nj = 2 flavours. 



the continuum limit. A similar effect has also been seen elsewhere, for example 
in the determination of roA MS for nj = flavours, where coarse lattices also 
show this characteristic flattening of the data. 

Finally, we compare these numbers with results from the QCD sum rule ap- 
proach. A review of results from this method is given in [S3], citing as a final 
result mf s (2GeV) = 99(28) MeV, while a recent five-loop calculation, gives 
m K s IS {2 GeV) = 105(6)(7) MeV. These numbers cover the lattice results in Fig.fTol 

In conclusion, although there is a spread of results, it would seem that the 
unquenched strange quark mass determined here is not lighter than the quenched 
strange quark mass and lies in the range of 100 - 130 MeV. 
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Appendix 



We collect here in Tables El El and El the numerical values of the partially 
quenched am ps and the bare AWI quark terms am^\ am^> and am q . These are 
defined in eq. ()39|) as ratios of certain correlation functions. (The operators are 
summed over spatial planes, the time derivatives being taken as in footnotes Eland 
El) The (bootstrap) errors for the ratios are given uniformly to two significant 
figures, with the overriding requirement that the result must also have a minimum 
of four significant figures. 

The second column in the tables gives the pion mass, defined in the standard 
way from the correlation function 

^psmeared^ psmeared t^t_0 ^ — m ps t _|_ ^ -m p3 (aN T — t) ^ (70) 

where the correlation function is evaluated in a configuration with sea quark 
mass arriq. (The normalisation is not important here so we can work with unim- 
proved operators to obtain the pseudoscalar masses.) The smearing used is Jacobi 
smearing (see e.g. |SZ|), with typical parameters k s = 0.21 and n s = 50. 

The third and fourth columns in the tables give the bare results for am^ and 
am'p, as defined in eq. (jSHJ)- 

The improvement coefficient ca has been determined non-perturbatively in 
|2"U] . We use the values obtained there of 



p 


ca 


5.20 


-0.0641(40) 


5.25 


-0.0565(40) 


5.29 


-0.0517(40) 


5.40 


-0.0420(40) 



(71) 



From the tables it can be seen that the inclusion of the improvement term (xca) 
to the quark mass gives a noticeable change in the final result. Also the error in 
ca has an effect. Although not a large difference to using simple error propaga- 
tion for the three quantities {arn^\ am^ and ca), to try to minimize the error 
propagation we have adopted the procedure of first finding the bootstrap error 
for am^ + CAorh^ (with fixed ca) and then including the independent error of 
ca (x fixed am' 1 ') by error propagation. This gives the results for am g shown in 
the fifth column of the tables. 



31 



Kg 


am ps 






arriq 








tv q — 


0.1342 






0.1334 


0.6581 


12) 


0.1122K 


19) 


0.22501(73) 


0.09778 


;i6) 


0.1338 


0.6224 


12) 


0.100101 


19) 


0.20041(70) 


0.08724 


16) 


0.1342 


0.5847 


12) 


0.0882K 


19) 


0.17647(67) 


0.07689 


16) 


0.1347 


0.5359 


12) 


0.07366( 


19) 


0.14744(63) 


0.06420 


16) 


0.1353 


0.4720 


13) 


0.05664( 


19) 


0.11383(60) 


0.04934 


;i6) 


0.1356 


0.4371 


14) 


0.04828( 


19) 


0.09743(59) 


0.04203 


16) 


0.1360 


0.3856 


16) 


0.03715( 


22) 


0.07576(60) 


0.03229 


19) 


0.1362 


0.3569 


17) 


0.03159( 


25) 


0.06515(64) 


0.02741 


[22) 








— 


0.1350 






0.1332 


0.5985 


11) 


0.10076( 


16) 


0.18529(60) 


0.08888 


;i3) 


0.1337 


0.5515 


11) 


0.08566( 


16) 


0.15670(58) 


0.07561 


;i3) 


0.1342 


0.5018 


11) 


0.07088( 


15) 


0.12923(55) 


0.06260 


;i3) 


0.1345 


0.4703 


12) 


0.06218( 


15) 


0.11327(54) 


0.05492 


;i3) 


0.1350 


0.4148 


13) 


0.04785( 


15) 


0.08755(54) 


0.04224 


13) 


U. looo 


0.3771 


iiS) 


0.03954( 


17) 


0.07246(54) 


0.03490 


;i4) 


0.1355 


0.3505 


19) 


0.03397( 


18) 


0.06256(57) 


0.02996 


;i6) 


0.1357 


0.3216 


;2o) 


0.02837( 


20) 


0.05293(57) 


0.02497 


;i7) 










0.1355 






0.1332 


0.5546 


10) 


0.0912K 


14) 


0.15890(54) 


0.08102 


;i2) 


0.1336 


0.5158 


ii) 


0.0791K 


14) 


0.13707(52) 


0.07032 


;i3) 


0.1340 


0.4751 


ii) 


0.06718( 


15) 


0.11600(51) 


0.05974 


;i3) 


0.1343 


0.4430 


12) 


0.05842( 


15) 


0.10069(50) 


0.05196 


;i3) 


0.1348 


0.3848 


14) 


0.04410( 


15) 


0.07610(51) 


0.03922 


;i4) 


0.1350 


0.3600 


liS) 


0.03847( 


15) 


0.06656(51) 


0.03421 


;i4) 


0.1353 


0.3200 




0.03014( 


16) 


0.05252(53) 


0.02677 


;i4) 


0.1355 


0.2907 




0.0245K 


15) 


0.04275(43) 


0.02177 


;is) 


0.1357 


0.2577 


|23) 


0.01909( 


18) 


0.03475(51) 


0.01687 


;i6) 



Table 7: The bare results for am ps , am q , am q and am q for (3 = 5.20. 
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Kg 


am ps 






am q 


Kg = 0.1346 


0.1337 


0.5794(15) 


0.09697(19) 


0.17225 


;87) 


0.08723(17) 


0.1340 


0.5514(15) 


0.08793(19) 


0.15554 


;86) 


0.07914(17) 


0.1346 


0.4932(10) 


0.07027(12) 


0.12399 


;47) 


0.06326(11) 


0.1349 


0.4612(17) 


0.06148(18) 


0.10782 


79) 


0.05539(17) 


0.1353 


0.4168(18) 


0.05005(18) 


0.08774 


76) 


0.04509(17) 


0.1355 


0.3932(18) 


0.04440(18) 


0.07794 


75) 


0.03999(17) 


0.1359 


0.3420(20) 


0.03317(19) 


0.05877 


;64) 


0.02985(18) 


0.1361 


0.3133(22) 


0.02750(21) 


0.04913 


;65) 


0.02472(19) 


«£' = 0.1352 


0.1337 


0.5419(11) 


0.08860(16) 


0.15186 


'60) 


0.08002(14) 


0.1341 


0.5027(12) 


0.07658(16) 


0.13073 


W) 


0.06919(14) 


0.1345 


0.4621(13) 


0.06474(16) 


0.11027 


;54) 


0.05851(14) 


0.1348 


0.4300(13) 


0.05599(16) 


0.09536 


;52) 


0.05060(14) 


0.1352 


0.3821(13) 


0.04432(12) 


0.07471 


;38) 


0.04010(10) 


0.1355 


0.3466(17) 


0.03593(17) 


0.06178 


jso) 


0.03244(16) 


0.1358 


0.3054(20) 


0.02740(19) 


0.04775 


;53) 


0.02470(18) 


0.1359 


0.2901(22) 


0.02452(21) 


0.04309 


;56) 


0.02209(20) 


Kq = 0.13575 


0.1336 


0.50970(72) 


0.084021(81) 


0.13257 


'36) 


0.076528(72) 


0.1339 


0.48011(72) 


0.074990(81) 


0.11733 


;35) 


0.068359(72) 


0.1343 


0.43883(73) 


0.063110(80) 


0.09772 


;33) 


0.057587(73) 


0.1346 


0.40619(74) 


0.054323(81) 


0.08354 


;3i) 


0.049601(74) 


0.1350 


0.35966(76) 


0.042776(83) 


0.06532 


;28) 


0.039085(76) 


0.1352 


0.33469(77) 


0.037071(84) 


0.05657 


I 24 ) 


0.033874(78) 


0.1355 


0.29421(81) 


0.028608(87) 


0.04360 


[22) 


0.026144(81) 


0.13575 


0.25556(55) 


0.021495(57) 


0.03291 


;i5) 


0.019635(52) 


0.1360 


0.2117(13) 


0.01456(11) 


0.02256 


;26) 


0.013281(98) 



Table 8: The bare results for am ps , am q , am q and am q for (3 = 5.25. 
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K q 


am ps 






am q 


K b q = 0.1340 


0.1340 


0.5767(11) 


0.09689(19) 


0.17170(60) 


0.08802(17) 


0.1344 


0.5392(15) 


0.08480(18) 


0.14841(82) 


0.07713(16) 


0.1349 


0.4901(16) 


0.07010(18) 


0.12204(81) 


0.06379(16) 


0.1352 


0.4589(17) 


0.06141(19) 


0.10669(80) 


0.05590(17) 


0.1355 


0.4255(17) 


0.05283(19) 


0.09169(80) 


0.04809(17) 


0.1357 


0.4024(20) 


0.04715(19) 


0.08189(80) 


0.04292(17) 


0.1359 


0.3781(21) 


0.04151(20) 


0.07226(80) 


0.03778(18) 


0.1362 


0.3384(23) 


0.03315(23) 


0.05830(81) 


0.03013(21) 


«£' = 0.1350 


0.1340 


0.52221(81) 


0.08542(11) 


0.14024(41) 


0.078169(94) 


0.1343 


0.49323(83) 


0.07643(11) 


0.12483(40) 


0.069978(96) 


0.1347 


0.45278(86) 


0.06460(11) 


0.10489(39) 


0.059180(99) 


0.1350 


0.42057(92) 


0.05584(11) 


0.09036(34) 


0.051171(95) 


0.1355 


0.3634(10) 


0.04146(11) 


0.06707(38) 


0.03799(10) 


0.1357 


0.3381(10) 


0.03577(12) 


0.05794(38) 


0.03277(11) 


0.1360 


0.2963(12) 


0.02722(12) 


0.04452(39) 


0.02492(11) 


0.1361 


0.2798(17) 


0.02427(15) 


0.03963(54) 


0.02222(13) 


= 0.1355 


0.1339 


0.49968(92) 


0.08260(11) 


0.12790(44) 


0.075988(94) 


0.1343 


0.46105(86) 


0.07062(10) 


0.10825(41) 


0.065025(93) 


0.1346 


0.43015(86) 


0.06175(10) 


0.09402(38) 


0.056888(92) 


0.1349 


0.39774(87) 


0.05297(10) 


0.08023(36) 


0.048825(92) 


0.1353 


0.35144(91) 


0.04142(10) 


0.06252(34) 


0.038192(93) 


0.1355 


0.32688(70) 


0.035783(77) 


0.05417(22) 


0.032983(71) 


0.1358 


0.2858(11) 


0.02723(11) 


0.04151(30) 


0.02508(10) 


0.1360 


0.2552(14) 


0.02159(13) 


0.03316(31) 


0.01988(12) 


0.1363 


0.2012(18) 


0.01304(14) 


0.02056(30) 


0.01197(13) 


K b q = 0.1359 


0.1339 


0.47757(65) 


0.078462(92) 


0.11620(29) 


0.072456(88) 


0.13425 


0.44247(66) 


0.067952(92) 


0.09946(28) 


0.062811(88) 


0.1346 


0.40540(70) 


0.057566(92) 


0.08333(27) 


0.053259(88) 


0.13505 


0.35469(70) 


0.044401(93) 


0.06350(25) 


0.041119(89) 


0.13531 


0.32287(73) 


0.036894(94) 


0.05251(24) 


0.034180(90) 


0.13562 


0.28151(78) 


0.028036(96) 


0.03982(23) 


0.025977(91) 


0.1359 


0.23924(87) 


0.020134(92) 


0.02886(21) 


0.018642(86) 


0.13617 


0.1899(12) 


0.01239(10) 


0.01837(25) 


0.011444(98) 



Table 9: The bare results for am ps , am q , am q and am q for (3 = 5.29. 
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Kq 


am ps 






am q 


= 0.1350 


0.1346 

0.1350 

0.1353 

0.1357 

0.13602 

0.1363 

0.13655 

0.1366 


0.44399(52) 
0.40301(43) 
0.37156(54) 
0.32541(63) 
0.28482(69) 
0.24504(77) 
0.20349(95) 
0.1934(13) 


0.071849(57) 
0.059913(50) 
0.051132(70) 
0.039519(68) 
0.030274(70) 
0.022220(72) 
0.014975(74) 
0.013555(84) 


0.10023(25) 
0.08250(18) 
0.06996(21) 
0.05359(21) 
0.04100(21) 
0.03027(20) 
0.02098(22) 
0.01950(39) 


0.067635(52) 
0.056444(47) 
0.048190(67) 
0.037266(64) 
0.028550(66) 
0.020947(68) 
0.014093(70) 
0.012735(79) 


«£' = 0.1356 


0.1346 

0.13494 

0.1353 

0.1356 

0.13591 

0.13618 

0.13643 

0.1365 


0.42009(66) 

0.38581(63) 

0.34617(72) 

0.31232(67) 

0.27210(77) 

0.23346(87) 

0.1921(10) 

0.1796(12) 


0.067720(51) 
0.057620(49) 
0.046987(49) 
0.038239(49) 
0.029197(52) 
0.021380(56) 
0.014081(68) 
0.01185(14) 


0.08973(29) 
0.07533(28) 
0.06068(26) 
0.04926(22) 
0.03725(24) 
0.02742(23) 
0.01873(22) 
0.01629(34) 


0.063947(47) 
0.054453(45) 
0.044436(46) 
0.036168(44) 
0.027631(50) 
0.020227(53) 
0.013293(65) 
0.01116(13) 


K b g = 0.1361 


0.1346 

0.13493 

0.13525 

0.13555 

0.13584 

0.1361 

0.13632 

0.1364 


0.40055(60) 

0.36621(63) 

0.33068(70) 

0.29521(76) 

0.25784(85) 

0.22081(72) 

0.1833(12) 

0.1668(19) 


0.064373(49) 
0.054572(47) 
0.045072(53) 
0.036298(51) 
0.027842(55) 
0.020335(47) 
0.013937(71) 
0.011552(88) 


0.08174(23) 
0.06820(22) 
0.05559(21) 
0.04425(20) 
0.03372(20) 
0.02455(16) 
0.01705(22) 
0.01457(22) 


0.060936(45) 
0.051704(44) 
0.042735(50) 
0.034437(49) 
0.026425(53) 
0.019303(44) 
0.013220(68) 
0.010940(84) 



Table 10: The bare results for am ps , am q , am q and am q for (3 = 5.40. 
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